Quasiparticles in the Vortex State of Dirty d-Wave Superconductors 
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We consider the problem of the vortex contribution to thermal properties of dirty d-wave super- 
conductors. In the clean limit, Volovik has argued that the main contribution to the density of 
states in a d-wave superconductor arises from extended quasiparticle states which may be treated 

semiclassically, giving rise to a specific heat contribution SC{H)'^,H 1 ^ 2 . We show that the ex- 
tended states continue to dominate the dirty limit, but lead to a HlogH behavior at the lowest 
fields, Hd <JJ« H C 2- This crossover may explain recent discrepancies in specicific heat measure- 
ments at low temperatures and fields in the cuprate superconductors. We further discuss the field 
dependence of transport properties within the same model. 



PACS Numbers: 74.25.Fy, 74.72.-h,74.25.Jb 

Introduction. As pointed out by Volovik, |ffl] d-wave 
superconductors are fundamentally different from their 
s-wave counterparts in the vortex state, in that at low 
temperatures and fields the extended quasiparticle states 
around the vortex play a much more important role. In 
the case of the single-particle density of states, he showed 
in fact that the contribution from these states over the 
entire vortex unit cell exceeds that of the bound states. 
It is therefore to be expected that the quasiparticles will 
also play an important role in transport properties of a 
d-wave system in the presence of an applied magnetic 
field. To study low-T transport, one needs a way of 
including scattering by defects at a simple level. We 



is the superfiuid velocity. The impurity-averaged electron 
propagator in particle-hole space is given by 



(Q - v s • k)r + A k Ti + £ k T 3 



(1) 



(S5-v..k)2-A£-g ' 

where the Tj are Pauli matrices. Due to symmetry of the 
d x 2_ y 2 order parameter A k = A cos 2c/> and the assumed 
particle-hole symmetry of the normal state only the fre- 
quency is renormalized uj = to — Y.q(lu). In the unitar- 
ity limit the self-energy is given by So = tqT/Gq where 
r = ni/irN is an impurity scattering rate depending 
on the concentration rii of point potential scatterers and 
the density of states at the Fermi level, Nq. For a simple 



recently studied the specific heat 6C(H) of d x i_ y i su- d x i_ y i superconductor the averaged integrated Green's 



perconductors in a field using the common <-matrix ap- 
proximation for strong impurity scattering, Q propos- 
ing a possible explanation for deviations of some mea- 
surements from the SC ~ VH behavior predicted for a 
clean d-wave superconductor by Volovik. We also showed 
how other consequences of the clean Dirac spectrum, 
including the l/u> divergence in the density of states 
5N(oj; H) H and the predicted scaling [|| of the spe- 
cific heat SC/iTH 1 / 2 ) = F C {H^ 2 /T), are modified by 
disorder. 

In this paper we review the results on thermodynamic 
properties and present a preliminary calculation of the 
low-field microwave conductivity. We neglect other pos- 
sible contributions to the conductivity arising from disor- 
dered vortex scattering, and assume that the Abrikosov 
lattice is pinned, ordered, and sufficiently dilute to not 
alter the quasiparticle bands significantly from their 1- 
vortex form. We point out that confirmation of the mag- 
netic field and disorder dependence would not only pro- 
vide evidence for the primacy of quasiparticle transport 
in d-wave superconductors, but also yield information on 
the phase shifts of the defects and on quasiparticle relax- 
ation times. 

Semiclassical treatment of extended states. In the 
semiclassical approximation the Doppler effect associated 
with supercurrents around a vortex lead to a shift of the 
quasiparticle energy u) — ► lu — v s -k, where v s = (h/2mr)9 



function reads Gq{uj, v s ) = — i(2/7r)K(Ao/(u> — v s • k)), 
where K is the complete elliptic integral of the first kind. 
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FIG. 1. V(0; H)/N for r/A = 0.1, 0.01, and 0.001 in uni- 
tarity limit (solid lines) and the clean limit (dotted line). Data 
from Fisher et al. |J (circles) ; Moler et al. |Q] (untwinned sam- 
ple, squares), assuming # c2 /a 2 =300T, 7„ = 15 mJ-mol-K 2 . 

Density of states at zero energy. To find N(0;H), 
we average the propagator over a vortex unit cell, 
N(0;H)/N = (-1mGo(u]v a )) B , where for any /(v s ) 
we define (f(v s )) H = A~ l {H) J ccll d 2 r f(v s ). The mag- 
netic field dependence enters only through the inter- 
vortex spacing R = £ (Tr/2) 1 / 2 a- 1 (H c2 /H) 1 / 2 . In 



1 



the clean limit we r eproduce the result of Volovik 
N(0;H)/N ~ ^/sJna{H/H c2 ) 1/2 . In the dirty limit, 
where the zero-energy quasiparticle scattering rate 70 is 
much larger than the average quasiparticle energy shift 
E H = a(iJ/ir c2 ) 1/2 A , we find 



N 8 7o \H c2 ' 



2a 2 I H 



(2) 



Numerical results for the density of states N(0, H) to- 
gether with experimental data of Moler et al. 0] and 
Fisher et al. are given in Fig. 1. The data Moler et 
al. are consistent with a slightly dirty d-wave supercon- 
ductor, while the Fisher data cannot be well fitted to the 
clean case, suggesting that their sample contains roughly 
an order of magnitude more defects. 

Density of states at finite frequency. The field depen- 
dent part of the density of states for v s kp <§C u> <C Ao 
is given by (-Im[G (w; v s ) - G q (uj)]) h f§ which for 
70, < lo yields 



C 



N{Q;H)\T 



N(uj;H) 



(4) 



T < max[7o,£; ff ] < A 
j Q ,E H < T< A 



A T 2 term characteristic of the pure d-wave system in 
zero field is present whenever both the impurity and mag- 
netic field scales are smaller than the temperature. Sub- 
stituting Eq. (3) into (4) leads to the scaling function 
of the specific heat 8C(H)/['y n Ta(H/H c2 ) 1 / 2 ] = F C [Y], 
where Y = a(H/H c2 ) 1/2 T c /T. The scaling function F c 
varies as F c ~ 31og2A F/(47rT c ) for 7 « 1 and as 
Fc — \J2pK for Y ^> 1. Scaling is expected for a given 
data set provided H, T are such that Eh and T are both 
larger than the impurity scale 70. A full numerical evalu- 
ation of (4) plotted in Fig. 3 shows that for the clean case 
(open symbols) , scaling is obtained over the full range of 
Y, whereas for the dirty system (filled symbols) scaling 
has broken down completely 
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'+(l + 2a; 2 ) 
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where x — \J2/tt(lu / Eh) and k„ the nodal direction. 

In Fig. 2 we plot the density of states versus frequency. 
In the intermediate range we recover the predicted l/u> 
divergence 6N(lj;H)/Nq ~ a 2 TrA H/(4:UjH c2 ), which in 
the clean limit 70 < Eh is cut off as SN(ui; H) ~ 
A(0; H)(l — 7rx/4)) and in the dirty case by the impurity 
scattering scale 70 > Eh ■ 
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FIG. 2. Density of states at a(H/H c2 ) 1/2 = 0.025 vs. freq. 
lu/Aq for r/T c =0., 0.0001,0. 001, 0.01 (solid lines). Interme- 
diate-frequency asymptotic result SN ~ a 2 TvAoH/(4:LoH c2 ) 
(dashed line). 



Scaling of specific heat. 
temperature is given by 



The specific heat at low- 



x > 1 
x < 1 



■ co 




d 



□ t=o.i r=o.ooi ■ t=o.i r=o.i 

o T=0.05 r=0.001 • T=0.05 r=o.i 
* T=0.025 r=0.001 * T=0.025 r=0.1 



X 

D 



X 

o 



CM 

d 




12 3 4 

a(H/H c2 ) 1 /2 Tc / T 

FIG. 3. Normalized vortex contribution to specific heat, 
5C e i(H)/[ lel Ta(H/H c2 ) 1/2 } vs. Y = a{H/H c2 ) 1/2 T c /T for 
fixed temperatures T and scattering rates T as shown; unit of 
energy T c0 . Asymptotic large- Y limit \J2fn (dashed line). 

Microwave conductivity. The microwave conductivity 
is given by [§ 



where 



ne 
m 



^sech 2 (^) • Sij(u,n) 



(5) 



Sij(u,Q) — Im / — kikj 
Ztt 



1 



1 



UJ+ — UJ_ 



1 



£0+ £,'0- 



(6) 



Here we have defined d>± — ui ~ ^o(lo ± «0 + ), £o± = 
± sgn(u))(u}|. — A 2 ) 1 / 2 , as well as analogous primed quan- 
tities uj± — uj±(u) — £1) and £ 0± = £,o±(u) — 0), with 
u> = to — v s • k„. Taking the limit f2 — > 0, T — > and 
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performing the (^-integration we find for the contribution 
to the conductivity from currents at a position r relative 
to the vortex center 



Su(u},t) 



1 



4 ^ ttA 

nodes 



Re<^ k 



-K(fc) - E(fc) 



(7) 



where K(fc) and E(/c) are the complete elliptic integrals of 
the first and second kind, respectively, with the argument 
k = A /(Aq — (S — Co) 2 ) 1 / 2 . Note the space dependence 
arises exclusively through the Doppler shift v s ■ k„ = 
(Tlvf/Zt) cos#. At low temperature and for H c \ <C H <C 
H C 2 the leading order is given by 



/ x o'oo 

nodes 



du -sech (— ) 




(8) 



where cr o = is the universal conductivity [|10| . 

In the gapless regime the self-energy reduces to — E ' ~ 
7o 3> which leads to the following leading order con- 
tribution of the magnetic field to the conductivity 



Sa(r) = (7 o 



To 



(9) 



In the clean limit 70 <C Eh the self-energy T, (uj = 0, r) 
is given by 



E o '(w = 0,r) 



Ao 



(10) 



that for small fields of order a few Tesla, vortices will be 
sufficiently dilute that scattering at low T will be impu- 
rity dominated. Orenstein et al. indeed observe a convex 
downward curvature in 5a(H), as well as a broad max- 
imum at around 3T. Further investigation of scattering 
of quasiparticles by vortices and of the effect of applied 
field on inelastic spin fluctuations is needed to extend our 
results to higher temperatures and fields. 

Note added: Results for N(0) in the dirty limit similar 
to ours were obtained independently by Barash et al. |l5| ] 
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which leads to 



Mr) = a„ - ( - E - 



Ao 

r 



(ii) 



By averaging the above results over a unit vortex cell 
we obtain the magnetic field contribution to the zero- 
frequency conductivity for a d-wave superconductor 
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70 < E H 



(12) 



These results should be valid provided the quasiparticle 
mean free path is in fact limited by impurities and not 
by fluctuations in the vortex lattice. The experiments of 
Orenstein et al. [Q at 150 GHz and a few Tesla are in 
precisely the correct regime to allow one to neglect ab- 
sorption into vortex oscillation modes, and to compare 
to our Q — > result. Furthermore, we expect even in 
the presence of elastic (pinned) vortex lattice disorder 
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